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Abstract 

We continue the study of U(\) vortices with cholesteric vacuum struc¬ 
ture. A new class of solutions is found which represent global vortices of the 
internal spin field. These spin vortices are characterized by a non-vanishing 
angular dependence at spatial infinity, or winding. We show that despite 
the topological Z 2 behavior of SO (3) windings, the topological charge of the 
spin vortices is of the Z type in the cholesteric. We find these solutions nu¬ 
merically and discuss the properties derived from their low energy effective 
field theory in 1 + 1 dimensions. 



1 Introduction 


This paper is a continuation of [ilj which investigates the low energy theory 
of U(l) vortices in systems with cholesteric vacuum structures. This work 
was inspired by considerations from supersymmetric solitons (see [2] for a 
review) as well as superconductivity |3] and superfluidity pE]. The kind 
of vortex solutions found in [1] are similar to the w vortices appearing in 
superfluid 3 He-B [lj. Indeed these w vortices, and those found in [T], develop 
a cholesteric spiral structure in the core. If a collection of such vortices are 
coherently orientated in a rotating cryostat the spiral structure appears 
in the bulk superfluid |4j. The main task of this paper is to report on 
a different class of solutions also present in this model. These solutions 
are characterized by the internal spin field carrying a topological charge of 
its own. Typically vortices with topological windings originating from an 
internal SO(3)j spin group are referred to as spin vortices. However, for 
reasons discussed below, the solutions we obtain are distinct from the typical 
spin vortex observed in say superfluid 3 He mm- Our plan is to find these 
solutions numerically and investigate their low energy excitations. 

The model, as originally introduced in [lj, has a symmetry structure 
given by 

G = 1/(1) gauge X SO{3)j X T (1) 

where the t/(l) gaug e group represents gauge rotations of a complex Higgs 
field, and SO(3)s+l = SO(3)j represents the spin orbit locked rotational 
symmetry. The group T represents the three translational symmetries of 
the model. The spin orbit locked SO(3)j appears from the larger group 
SO(3)s X SO(3)l due to an additional term in the Largrangian coupling 
the spin index to the spatial gradiant. This term is parity violating and 
leads to spontaneous breaking of translational symmetry. The additional 
term thus gives rise to a cholesteric spiral structure in the vacuum under 
certain conditions of the parameters in the model [I]. 

The spontaneous breaking of the U( l) gaug e group allows for the ap- 
pearence of the Abrikosov-Nielsen-Olsen [7] flux tube familiar to models 
of type II superconductors [3] and Yang-Mills theories [2J. With no other 
winding of the subgroups of G these vortices represent the mass vortices 
that were studied in [T]. A typical spin vortex on the other hand arises from 
the non-trival topology of the SO(3)j group: 

vn(50(3)j) = Z 2 . (2) 

The spin vortices can appear regardless of the cholesteric behavior of the 
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vacuum and only require the existence of a broken SO( 3) symmetry. How¬ 
ever, the vortices we consider in this paper are distinct from the typical 
spin vortex in that they are only allowed due to the broken translational 
symmetry of the vacuum. In particular, the vacuum of our model preserves 
a U(l)j, global subgroup of G which is generated by a linear combination 
of the translational and rotational generators about an axis oriented along 
the cholesteric wave vector k: 

G -» 1/(1)J/ x T X) y, where 17(1) j' C U(l)j z x T z . (3) 

T xy is the translational group along the direction perpendicular to k, whilst 
T z is that along k. Here the subgroup 1/(1) j/ is denoted with its generator 
J' z which is a linear combination of the rotational J~ and translational p z 
generators: 

J'=J z -j, k = \k\. (4) 

As we will discuss below this equivalence of rotations and translations in 
the vacuum leads to a global U( 1) degeneracy space of vacua in addition to 
the U (1)g aU ge “degeneracy” appearing from the Higgs phase. Vortices may 
appear with either gauge or global U( 1) topological windings or both: 

7r l(//(l)gauge X 1/(1)) = Z + Z. (5) 

In this paper we will focus our attention on the vortices with both types 
of charges. We hasten to point out that the vortices in this model are global. 
Thus, for cholesteric vacua (vacuum I as it is defined below) these vortices 
will have a logarithmically divergent energy functional and thus require an 
infrared cutoff. This is not a problem for systems in a finite volume such as 
superfluid 3 He. However, as a consequence of this finite volume, the Kelvin 
modes appearing on the vortex will acquire a mass gap proportional to the 
inverse cutoff length R. These effects were considered recently in [8], and 
we will apply similar arguments when we discuss the gapless excitations in 
section 5 below. 

We will begin in sections 2 and 3 with a detailed review of the model 
and the vacuum structure following the procedure discussed in [T|. We will 
explore briefly the Goldstones appearing in the bulk cholesteric vacuum. We 
will discuss their relation to the spin vortices appearing in the cholesteric 
bulk. In section 4 we will present the equations of motion for the spin vortices 
and present their numerical solutions for a representative set of parameters. 
In section 5 we will discuss the low energy theory of the moduli appearing 
on the vortices. We will present an analysis of the Kelvin modes for vortices 
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in a particular vacuum, and provide additional topological discussions of the 
orientational moduli. For the case of a cholesteric vacuum we will discuss 
the finite volume effects on the mass terms for Kelvin modes. Finally, we 
provide concluding remarks in section 6. 


2 The system 

The starting point is the system used in [1J. This model is an extension of 
the model originally suggested in [9], inspired by Witten’s superconducting 
string model fOj. and further studied in m- Its Lagrangian has the form 


C, — Cq + jC y T C-e 


£° = + \D^\ 2 - A(H 2 - v 2 ) 2 , 


C-e 


dtXidtXi - diXjdiXj 
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(-^ 2 + \(j>\ 2 )XiXi + PUiXi ) 2 , 


V^ijkXi^jXk j 


( 6 ) 


where the subscript i runs over i = 1, 2, 3. The field Xi can be viewed as a 
spin field. The covariant derivative is defined in a standard way 


D^ = d^~ iA^ . (7) 

The last term in Q is the parity violating twist term with a single derivative. 
This term is responsible for the cholesteric behavior of the vacuum as we 
will show below. 

Let us consider the system with r/ = 0 first. Then, under an appropriate 
choice of parameters the charged field </> condenses in the ground state, 


I vac — ^ i 


( 8 ) 


and, if ji < v, the field % is not excited in the vacuum 

te)vac = 0- (9) 

The system has full translational invariance T and orbital rotational sym¬ 
metry SO(3)l as well as an internal rotational SO(3)s symmetry of the x.i 
spin field. As is well known, the model supports the Abrikosov flux tube [?]. 
Inside the tube, where the Higgs field vanishes, the spin field Xi is excited 


3 




giving rise to gapless (or quasigapless) excitations of non-Abelian type, lo¬ 
calized on the vortex. The full solution is known as a non-Abelian string. 

Including C £ , explicitly breaks the orbital rotational part of the Lorentz 
symmetry implying a spin-orbit locked symmetry of the full Lagrangian, 

SO(3)l x SO(3) s ->• SO(3) s+ l = SO(3)j. (10) 

The energy density derived from the Lagrangian ([h]) is 

E J 2 /■'./ /' ^ + I A^l 2 + xm 2 - V 2 ) 2 + diXjdiXj 

+ V£ijkXidjXk + 7 [(-E 2 + \(t>?)XiXi + PiXiXi?] ■ (H) 

Note that the kinetic terms of xi including the linear in derivatives rj term 
can be recognized as the Frank-Oseen free energy density of an isotropic 
chiral nematic liquid crystal |12j . 

For later convenience we rescale all couplings in the model to make them 
dimensionless. Hence, we take 


z = , 


( 12 ) 


to denote the distance in the longitudinal direction aligned with the string 
axis. We will also set the dimensionless radial distance 

p = m^x 2 + y 2 , (13) 

where the mass of the elementary excitations of the charged field 4> is 

= 4Au 2 . (14) 

Additional dimensionless parameters are given as 
7 (c- 1) 


b = 


4Ac 


v 2 e 2 

c = A ’ a= 2X , f ) = T l/ r n<t>' 


(15) 


We limit ourselves to the semi-classical approximation and assume all cou¬ 
plings in the model are small. 

The field Xi being represented in Cartesian coordinates takes the form 


Xi = 






[xx(x,y,z), Xy(x,y,z), Xz(x,y,z) }. (16) 
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The static classical equations of motion are derived by extremization of 
energy in a general coordinate system they read 

dV 

V,: {V^99 l ^jXk) ~ VZkjrfjXi ~ 7 ^ = 0, 

A {g^V—gDjcj)) + (2A (|</>| 2 - v 2 ) + 'yg lJ XiXj) <\> = 0, 

di (V=99 ji 9 kl F jk ) - ie 2 V=~g (</> <*D l cf> - cj>D l ct>*) = 0 , (17) 

where 

f)y 

= 'T z gi {{-g 2 + \<t>\ 2 ) + 2/3 g lJ XiXj) Xk, (18) 

and 

ViXj = diXj ~ r \jXk (19) 

is the standard curved space covariant derivative with T the Christoffel sym¬ 
bol. 

3 Ground state 

In this section we give a summary of the results obtained in [1J concern¬ 
ing the ground state structure of the model. As shown there this model 
has a rich ground state structure depending on the value of and the re¬ 
lation between the parameters. Given the presence of the first derivative 
term one expects that in the ground state translational invariance will be 
spontaneously broken. If one assumes that this breaking is aligned in the z 
direction, then a cholesteric structure appears in the ground state with the 
X held rotating in the (x, y ) plane 


Xi = Xo ei(z), 


VW 
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Xi = Xo €i(z) , 


e(z) = |cosfcz, sin kz, oj 


( 20 ) 


where k is the z component of the wave vector k. xo is the constant ampli¬ 
tude of the x field in the vacuum to be determined below in (22). Then, in 
the ground state we may take 


<t> = , A = o, 


( 21 ) 
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where </>o is a constant. The equations of motion © imply 

Xo (- k 2 + r)k - 7((-/i 2 + <$j) + 2/3x<j)) = 0 > 

«/»o(2A(^-u 2 )+7x 2 ) =0. (22) 

Minimizing over the wavevector k one obtains 

k = rj/2. (23) 

Therefore, only when ij = 0 we recover translational invariance (however see 
( p0| ) below). As per |I) we always assume 

c>l, b > 0, j3 > 0 and a > 0, (24) 


and we use here the same dimensionless field definitions 


Xo 

<P 


rfH 

o 

yZ 


V2J3 X ° 

V(j). 


(25) 


There are two parameter ranges which characterize the vacuum struc¬ 
ture. We will report here the main results and refer the reader to [I] for 
calculational details. For both cases the lowest energy vacuum solution is 
dictated by the value of fj. This is important in order to understand the role 
of the spin-vortex solutions and their winding at spatial infinity. Consider 
first the case where 

/3(c - 1) > be. (26) 


Then when 


46 < r) 2 < 4/3 ( 1 — 


the vacuum solution of (22) is 




(27) 


HY = 


i - 


0(c- 1) 


rj 

4/3 


1 - 


be 

/3(c—1) 


(Xo) 7 = (c - 1) 


nl _ i 

46 1 


1 - 


be 

/3(c-l) 


(c — 1 ) 2 (fj 2 — 4b) 2 
/3(c — 1) — be 


(28) 
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In this vacuum both cpo 7 = 0 and xo 7 ^ 0 and therefore both local 17(1) gauge 
and z—translational symmetries are broken. This vacuum allows for the 
existence of local ANO vortices. If instead 


t < Ci tl = 4 & (29) 

then the vacuum solution is 

(¥>o)" = l, {xo) H = 0, EH C = 0, (30) 

Note that, since (xo) H = 0, this vacuum preserves translational invariance. 
Finally, when 

??crit 2 = 4 /3 ( 1_ ^(c-i)) ’ (31 ) 

the vacuum solution is given by 

M m = o, (xo) m = 1 + (c- 1)^, 

-in = 16bc(c - l)/3 - ((c - l)f ) 2 + 8 fe ) 2 
mc 646c(c - l)/3 

In this vacuum translational invariance is broken but local 17(1) gauge in¬ 
variance is preserved, hence this vacuum does not support ANO vortices. 
The energy of the ground states as functions of r] are shown in Figure 1. 

In the opposite case where 



/3(c — 1) < be, 


(33) 


the minimizing vacuum is given by (30) when 


~2 . ~2 
V ^ %rit3 


y/(c— 1 )bc/3 — b 

(^i) ' 


(34) 


When f) crosses the critical value r)^ rit3 , there is a first order transition from 
vacuum (30) to vacuum (32). This pattern of ground states is depicted in 
Figure 2. 

For this paper we will focus our attention primarily on the vacuum I and 
vacuum II cases where the gauge I7(l) gauge group is broken by a non-zero <f> 
field. For vacuum II none of the global symmetries SO(3)j or T are broken 
in the ground state. For ANO vortices in this case, cholesteric behavior of 
the Xi field will only appear in the vortex core. On the other hand vacuum / 
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has an interesting symmetry structure appearing from a linear combination 
of the rotational J z and translational p z generators: 


j> -J - Pi 


(35) 


From (20) for non-zero xo it is clear that J' z annihilates the vacuum, and 
thus the vacuum I ground state retains a £7(1) j/ symmetry in addition to 


the subgroup T xy of translations perpendicular to k = kz. We may thus 
write the symmetry group of the vacuum / 


Hi = U{l)j' z x T xy . (36) 

The equivalence of translations and rotations about the z axis on the 
vacuum state implies an additional structure on the degeneracy space of 
vacua in vacuum I. Consider the coset space 

G/H, = £7(1) gauge X Sj ± X U(l)j»/Z2, (37) 

where 

J” = J Z + ^ (38) 

is a generator of the £7(l)j// group orthogonal to J' z . Here Sj ± is the two 
sphere of directions for the wave vector k. The group £7(1) generates 
non-trivial rotations of the vacuum about the wave vector k. Physically the 
group represents degeneracy associated with rotations of x about k. 

Alternatively, due to the equivalence of translations and rotations, and the 
periodicity of the 2 : translations of the vacuum, the £7(1) j// may be viewed 
as a compactification of T z . 

Naively the topology of S 2 x U( 1)/Z2 supports only a Z 2 winding. How¬ 
ever, it was shown in P33 that for low energy excitations the generators J x 
and J y acquire a mass gap through an effective Higgs mechanism. For the 
vacuum state this implies that only a single Goldstone mode associated with 
phase rotations of the vacuum appear, instead of the naively expected three 
Goldstones from the three broken generators. Effectively we may write 


G/H, ->• H(l) gauge x U(l)jn. (39) 

The effective U(1)j>j degeneracy of vacuum / implies an additional topolog¬ 
ical charge associated with the integer winding of a closed contour about 
the vortex axis. We will explore numerical solutions with this winding in 
the next section. 
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Figure 1: b = 1, c = 1.25, f3 = 8. Vacuum energy dependence on fj for 
/3(c — 1) > be, the solid line corresponds to vacuum II, the dashed line to 
vacuum I and the dotted line to vacuum III. 


E 



Figure 2: 6=1, c = 1.25, {3 = 4. Vacuum energy dependence on fj for 
(3(c — 1) < be, the solid line corresponds to vacuum II and the dotted line 
to vacuum III. 
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4 Winding solutions 


In pp the full set of eight coupled second order three dimensional equations 
was solved with a pseudo-spectral numerical procedure, additional details 
of this procedure are provided below (see [33] for a detailed review of this 
method). The solution found, in which the spin field % asymptotes to the 
vacuum at spatial infinity, describes a local 17(1) Abrikosov vortex with a 
cholesteric spiral inside its core. These vortices resemble the w vortices 
appearing in superfluid 3 He-B [3]. 

The degeneracy space 17(1) j" allows for another class of solutions similar 
to the “spin-mass-vortices” (see ED with a global 17(1) winding, charac¬ 
terized by a dependence of the x.% field on the polar angle 9 at infinity. The 
8 dependence of Xi a t infinity is equivalent to translations of the vacuum 
state in the k direction. The compactified translational group allows for the 
topological integer winding about a closed spatial contour. 

As a consequence of the global nature of these vortices the profile func¬ 
tions at large p are IR divergent and require a regularization. Nonetheless, 
these solutions are important in the context of superfluid 3 He [5] . We are 
interested in solutions with spontaneous breaking of gauge symmetry, hence 
we only consider vacua I and II below. 

The geometry of our system is shown in Figure 3. In order to find the 
17(1) winding solutions we switch to axially symmetric coordinates and 
use the ansatz 


Xi 

Af 

<t> 


^=(Xr(p, z),xe(p, z), Xzi.Pi z )), 
A, = 0 , A e = 1 — f(p ), 
v<p(p) e ld ■ 


(40) 

(41) 
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Then extremization of the energy density gives the following equations, 


1 


P 


-d P ( pdpXz) + d'ixz ~ fj ( dpxe + — ) -^r(c^ + X 2 - 1 )Xz = 0, 


P 


c — 1 


'• ''' b -(cip 2 + x 2 ~ l)xe = 0 , 


-dp ( pdpXB ) + d z xe - -2 - »7(3*Xr - <9pXz) - 

p p A 

-dp ( pdpXr ) + d\xr - - 


c — 1 


P" 


c — 1 


(c<P + x - l)Xr = 0, 


- p dp{pd p y) + d 2 z ip - 




X 


<P = 0 , 


d 2 p f ~ ~ p d p f + d 2 f - cup 2 f = 0. 


(42) 


We wish to solve the equations (42) with the following boundary conditions 


<p(0) 

= 0, 

<p(oo) = <p 0, 


/( 0) 

= 1, 

/(oo) = 0, 


x'z( 0,z) 

= 0, 

O 

II 


Xr(0,z) 

= 0, 

Xr(oO,z) = X0 COS ( 

T). 

Xe{0,z) 

= 0, 

Xe(oo,z) = xo sin 1 

T). 


(43) 

where xo and (p are the same as the previous section and we demand con¬ 
tinuity in the z direction. In terms of the cartesian vector this boundary 
condition takes the form 


Xi = Xo{ cos (kz + 9), sin (kz + 9), 0} 


(44) 


which is not to be confused with the 17(1) vortex solution of (lj. In fact, it 
is only because of this subtle difference between asymptotic conditions that 
the full set of three dimensional equations © becomes two dimensional. 
This greatly simplifies the numerical procedure. 


4.1 Outline of numerical procedure 


The reduction of the spatial dependence of the coupled system (42) from 
three to two coordinates greatly simplifies the numerical analysis by reducing 
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Figure 3: Geometric set-up of the problem. The vortex axis, in the z di¬ 
rection, is parallel to the wave-vector k which is normal to the cholesteric 
planes. 


the computing power needed. However, for precise numerical results, two 
dimensional finite difference methods are still impractical when efficiency 
and precision are the goals. To achieve these goals we chose instead to use 
a pseudospectral method on a collocated grid. In this section we will give a 
brief outline of our procedure. A detailed and very readable description of 
these methods can be found in 

For the z dependence of the functions defined above, we choose to ex¬ 
pand in Fourier modes and evaluate the functions at the discrete points 
z m = 2irm/kM, where M is the total number of Fourier modes used in the 
decomposition. To decompose the r dependence in the functions above, we 
first compactify the radial direction using the transformation: 


Lx 

\/l — x 2 ’ 


x G [0,1]. 


(45) 


The parameter L determines the distribution of grid points. For this proce¬ 
dure we selected L = 2 to achieve the desired density of grid points near the 
origin. The radial dependence of the functions is then expanded in rational 
Chebyshev polynomials T n (x{r)) on a collocated Radau grid: 


x p 


= sin 


ir(2p — 1 ) 

AN — 2 ’ 


n G [1 ,N], 


(46) 
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Figure 4: The graph shows the solutions for <j>{r) (solid) and Aq(t) (dashed) 
at z = 0, for r) = 2.5. These solutions have negligable z dependence. 


where N is the total number of Chebyshev polynomials used in the decom¬ 
position. The Radau grid allows for boundary conditions to be imposed 
at spatial infinity with no conditions required at the origin. For functions 
required to vanish at the origin we expand in odd Chebyshevs. 


The pseudospectral decomposition allows the equations (42) to be writ¬ 


ten as a set of algebraic equations to be solved for the coefficients in the 
expansion. For the five functions to be solved for above, this amounts to 
solving a coupled system of 5 x Nx M polynomial equations. To approximate 
a solution the system is linearized and solved recursively. Each iteration of 
the procedure is repeated until convergence is achieved. 


4.2 Vacuum I 

Let us first consider solutions in which xo 7^ 0- I n this case the vacuum 


of the system is vacuum I, given by (28). Clearly in this case the solutions 


with boundary conditions (43) don’t approach the vacuum at infinity and 


are logarithmically divergent in the IR. These solutions require an IR reg¬ 
ularization. We restrict our solutions to a finite radial distance L. The 
winding of these spin-vortices ensures their topological stability. It should 
be pointed out that near the vortex core a non-zero value of Xz appears due 
to the variation of xe i n the radial direction. 

For parameter values a = 1, 6=1, c = 1.25, /3 = 8, and r] = 2.5 we 
obtain the solutions shown in Figures 4-6. 
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Figure 6: Shown are the vector fields x a t various values of the coordinate 
2 ; G [0,37r/4fc] for r] = 2.5. We note that as the vortex core is approached a 
non-zero value of Xz appears. This is expected from the equations of motion 
when x.O has an r-dependence near the core. 
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Figure 7: The graph shows the solutions for (f>{r) (solid) and Ag(r) (dashed) 
at z = 0, for r) = 2.2. These solutions have negligable z dependence. 


4.3 Vacuum II 


When xo = 0 the vacuum has the structure of vacuum II (30). In this 


case, the asymptotics of the solutions (43) have no cholesteric structure 
and thus we expect finite energy solutions. However, the winding nature 
of these solutions is lost. These solutions lack the topological stability of 
the previous case in vacuum I. They are instead metastable solutions whose 
topological stability is only supported to finite radial distance from the core. 
Presumably, these solutions will decay to the zero winding state through 
quantum tunneling via bubble nucleation. We leave an investigation of the 
tunneling rate to another project. Figures 7-9 illustrate these solutions for 
a = 1, 6 = 1, c = 1.25, /3 = 8, and r/ = 2.2. 


5 Low energy theory of vortex moduli 

To complete the analysis of the spin vortices in vacuums I and II, we discuss 
the low energy dynamics of their excitations localized on the vortex core. 
In [lj the low energy expansion was done at least for small r] for both the 
translational (Kelvin) and orientational moduli appearing in vacuum I and 
II. Here we will do an analogous treatment for the case of vacuum I spin 
vortices only. The moduli appearing in vacuum II are difficult to analyze 
due to the coupling of all translational and orientational modes, and give 
little additional physical information. We will content ourselves with a brief 
topological analysis of the vacuum II case, leaving the details for future 
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Figure 9: Shown are the vector fields x at various values of the coordinate 
2 G [0,37r/4fc] for the critical point rj = 2. We note that as the vortex core 
is approached a non-zero value of appears. This is expected from the 
equations of motion when %e has an r-dependence near the core. 
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work. 

It is best to consider topological arguments at first to provide a guide 
in the search for the gapless excitations appearing on the vortex core. We 
begin with vacuum I where the non-zero value of the spin field Xi breaks 
the translation symmetry in the z direction. The global symmetry group of 
the lagrangian ([b]) is given by 

^global SO(3)j x T (47) 

where T represents the three dimensional translational symmetry. In vac¬ 
uum I the ground state is invariant to the symmetry group 

Hr = *7(1)4 X Hxy C Gglobal- (48) 

However, in the asymptotic limit r —> oo the vortex posseses complete axial 
symmetry 

H' = *7(1)4 x T xy . (49) 

Gapless excitations of the vortex must appear from generators of H\ that 
are broken in the core. All other modes appearing from the generators of 
Ggiobai would have infinite tension in the kinetic terms. They would not be 
normalizable. 

In the previous section we found vortex solutions by assuming rotational 
symmetry about the z-axis. Thus, since the vortex solutions automatically 
break the translations in the x and y directions, our vortex solutions obey 
the symmetry group 

H'l, vortex = U(l)j z . (50) 

Thus, the degernacy space of the vortex solutions in vacuum I follow from 

tfj/tf/,vortex = (*7(1)4 X T xy )/U(l)j z = T xy . (51) 

We should thus find two moduli fields £ x>y (t,z) associated with the broken 
translational symmetries by the vortex axis. These moduli will however 
acquire a mass gap m ~ 1/R due to the IR cutoff imposed at large distances 
R from the vortex core. The details will be shown below. 

For the case of vacuum II the cholesteric structure does not appear in 
the vacuum and thus the full symmetry group of vacuum II is given by: 

Hji = Ggi oba i = SO(3)j x T. (52) 

As in the case of vacuum / the vortex preserves the subgroup SO{2)j z thus 
the degeneracy space is: 

Hn/Hn ,vortex = (SO(3)j X T)/U(l)j z = Sj ± x T. (53) 
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In this case the cholesteric structure only appears in the vortex core and thus 
an additional longitudinal translational mode appears on the vortex. The 
expansion of the low energy Lagrangian in terms of the gapless modes in¬ 
troduces additional interactions between the translational and orientational 
modes. We have omitted this topic for sake of clarity and space. 

We proceed with the analysis by expanding the Lagrangian by the small 
adiabatic translations £ x ,y(t,z) of the vortex solution 

<t>{x) - £(t, z)) 

An(x) ->• Afj,(x — £(t, z)) 

x{x) -> x(*- £(*>*))■ (54) 

Here £(t,z) only oscillates along the directions x and y perpendicular to 
the vortex axis. Preceding with the expansion of we arrive at the low 
energy theory 

5 2 E = S ab d z UU b + Ae ab tad z £ b , (55) 

where a and b are indices for the x and y directions. The symbol e ab = —e ba 
is the two dimensional antisymmetric symbol. The parameters S ab and A 
are given by 


S ab = / d 2 x±d a Xid b Xi, 


1 


A — r)^ab I d X± {^O z d a \id b \i + !j( run () (I \ n () b X/m ) • 


(56) 


Here m and n run over the x and y indices. We point out that S and A 
are in general functions of the z coordinate, owing to the dependence of the 
spin held x l ° n z - 


The expansion in (55) considered only the kinetic terms in £ a (t,z). We 


hasten to point out that this expansion is only valid up to the infrared 
regularization of our vortex solution far from the core. However, due to the 
global nature of the spin vortex, the hard cutoff of the vortex solution in 
the infrared leads to a non-local boundary contribution to the low energy 
theory. The infrared cutoff explicitly breaks the translational symmetry in 
the x and y directions and thus our low energy theory should include mass 
terms proportional to the inverse cutoff m ~ 1/R: 


<5 2 Ar = S ab £, a tib - S ab d z i a dzih ~ Ae ab ^ a d z £, b - rn 2 £ a £a- 


(57) 
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The details of the calculation of m 2 are given in the appendix. For a cylin¬ 
drical boundary we find that 


m 


2 


—7r 


Xo 

R 2 ' 


(58) 


and hence the mass term is negative implying an instability of the vortex 
solution. We should point out however, that this value for m 2 is dependent 
on the shape of the boundary. For non-rotationally symmetric boundaries, 
the mass term will in fact be tensorial: 


m 


m ab- 


(59) 


We consider only the rotationally symmetric boundary. 

For relativistic lagrangians this instability means the vortex wants to es¬ 
cape the finite volume. For non-relativistic cases, the negative mass implies 
a threshold for the dispersion relations where the Kelvin modes change their 
direction of propagation |8j. 


For the equation (57) the more appropriate context to consider is the 


relativistic case - although this is not strictly a relativistic system because the 


third term in (57) is not Lorentz invariant. There are two important limits 


in this system characterised by the wavelength of the modes A compared 
to that of the cholesteric structure fc . In the limit in which A > 1 /k the 
asymptotic values of the x % Held dominates the integral appearing in (56), the 


^ dependence of S a b disappears and we have that A = 0 identically. The long 
wavelength modes effectively ignore the “high” k structure of the vortex. As 
A ~ 1/k the wavelengths become comparable and we may no longer ignore 
the z dependence of S a b . Clearly, d ^ 0 in this limit. This region becomes 
difficult to tackle analytically and we won’t consider it further. We may 
however approximate an intermediate region in which A > 1/k where the z 


the isotropicity of the matrix S ( S xx = S yy and S xy = S yx 


dependence of S a b and A is small but 4 / 0. Then, varying (|57|) and using 
the isotrc 
equation 


Sac (fa - fa) + ~ m^acfa = 0. 

Taking the ansatz 

fa = C^+^fa, 

where f a is a constant vector, we may reduce this to 


) we obtain the 
(60) 


(61) 


M a bS,b = 0 , 


(62) 
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with 

M ab = (s ab (u : 2 - k 2 ) + Aie ab k - m 2 d ab ^ , (63) 

for the vector components of £ a . Finding the corresponding inodes is there¬ 
fore equivalent to finding the zero eigenvalues of the matrix M ab . The eigen¬ 
values are 


C*4 = (w 2 - k 2 )Syy - m 2 ± \J ( Ak ) 2 + ( {k 2 - u 2 )S yx ) , (64) 

setting these to zero gives four solutions for u, 


w? = 


= 


k 2 i s yx - s yy) ~ m2S yy + \j{ A kS yy ) 2 + S 2 X JrrA - ( Ak)‘ 


02 _ C2 

u yx u yy 


k 2 (Syx - Sly ) - m 2 Syy - y (AkSyy ) 2 + S 2 X (m 4 - ( Ak ) 

02 _ 02 
u yx u yy 


The (un-normalized) eigenvectors are 

t ( iyj(Ak) 2 + ((k 2 - u 2 )S yx ) 2 
^ ~~ 1 ± Ak - i(k 2 - uj 2 )S yx 

which define the diagonal basis as 

Ak — ( k 2 — oj 2 )S yx 


£diag I 


{£x £,y)i £x T £y 


(65) 


( 66 ) 


(67) 


( 68 ) 


to 


'-(Ak) 2 + (k 2 -u 2 ) 2 S* x 

The limit A> 1 /k (this k should not be confused with k) corresponds 

(69) 


2 ’ 7? 


Sab VTXo log — 5 ab , and A 0. 
ro 


where ?’o represents the coherence length of the vortex core. In this limit 
the system is Lorentz invariant and we find the dispersion relation: 


LO 


= ±\ k 2 — A 2 , where A 2 = 


1 


R 2 log R/ro' 


(70) 
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Here 2ir /A = A c is the critical wavelength of stability. Modes with wave¬ 
lengths longer than A c have purely imaginary u and are thus tachyonic 0 
For the specific case k = 0 the mode represents a complete transverse trans¬ 
lation of the vortex solution. Thus the vortex tends to escape from the finite 
volume. On the other hand wavelengths shorter than A c are real and thus 
represent stable solutions for the vortex [8j. 

In the proximicity of the long wavelength limit, where we may still con¬ 
sider equation (60) as a good approximation for the low energy modes, higher 
order terms in the S a b and A calculation need to be included. Even though 
we don’t perform this computation explicitly we note the main important 
feature of having a non-zero A: it introduces terms in the dispersion rela¬ 
tion of the form oj oc which is characteristic of surface waves in deep 
water. Finally, leaving the long-wavelength limit entirely equation (60) can 
no longer be trusted and one is forced to resort to numerics. We leave this 
treatment to a future publication. 


6 Conclusions 

The goal of this paper was to consider an extension of the analysis presented 
in [I] where mass vortices with a t/(l) gauge charge were considered in sys¬ 
tems with cholesteric vacuum structure. Such vortices play a role in many 
contexts including supersymmetric solitons [2], type II superconductors 0, 
and superfluid 1 * 3 He pi]. In the current paper, we extended the analysis to in¬ 
clude a different type of vortex, which supported a global topological 
charge in addition to the U( l) ga uge charge. Spin vortices of a similar descrip¬ 
tion were first observed in superfluid 3 He §3]. However, for the cholesteric 
vacuum state we have discussed above, the equivalence of rotations and 
translations along the z direction in the vacuum leads to a reduction of the 
degeneracy space SO(3)j —> U(\)jn. 

Due to the cholesteric behavior of the vacuum I solution, the U(1)j>j 
vortices are stable to deformation E3j. The low energy effective held the¬ 
ory of these vortices thus gives rise to a 1 + 1 dimensional held theory of 
translational excitations. A similar result is true in vacuum II, however the 
solutions here are free to tunnel to solutions with no U(1 )jn charge since 
the asymptotic form of the vacuum has no cholesteric structure. We have 
left a calculation of the tunneling process to future work. 

1 The “trivial” solution in which the imaginary part of uj has a sign corresponding to 

decaying modes represents a short-lived decaying mode in which the vortex core remains 

in its original position. 
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We have focused our attention on the classical description of the 1 + 1 
dimensional effective field theory of gapless excitations. Specifically, we’ve 
calculated the low energy Lagrangan for Kelvin modes of vortices in vac¬ 
uum I in the long-wavelength limit. For a relativistic dispersion relation a 
classical analysis is all that is necessary to count Goldstone modes appear¬ 
ing on the vortex. However, for the non-relativistic despersion relation, the 
quantization of modes becomes more complicated. Some general results for 
non-Lorentz invariant systems are discussed in m, and to systems with 
broken spacetime symmetries in fTS], however we have omitted a detailed 
analysis here. We have however explored the effects of a finite volume on the 
gapless excitations appearing on the vortex as was first pointed out in fRj. 
In particular, the global behavior of the vortex solution in vacuum I lifts 
the translational degeneracy, and mass gaps are generated for the Kelvin 
excitations. 

For vacuum II additional orientational modes appear on the vortex, 
however the details of these modes have added complications and we have 
left their analysis to future projects. 
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Appendix 

In this appendix we wish to show the details of the calculation of the bound¬ 
ary terms leading to the mass term appearing for the translational moduli 


in (57). We will first show that the linear term in vanishes, confirming 


that no instabilities appear on the vortex. Following this calculation we will 
expand to second order in £* and determine the mass term for the particular 
vortex configuration. 

We begin with a calculation of the linear term for the expansion pa¬ 
rameter £, which we may take to be constant. The t and z dependence 
is accounted for in the kinetic terms shown in (55). To first order in 
SXi = -iadaXi 

dv 

5E = 2di(6xjdiXj) + ri£ijk(SxidjXk + XidjSxk) + dxizz— 

OXi 

= 2 di (SxjdiXj ~ ijkXj^Xk ) 

+2 Sxi (~d 2 Xi + V e ijkdjXk + 

= — 2 f a di (daXjdiXj T 2 £ ijk^aXjXk'j i (71) 


where the second line in the second equality vanished by the equations of 
motion for Xi- To obtain the contribution of the linear term to the low 
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energy theory, we integrate over the perpendicular directions x and y: 


SE 


2£a J d~x±di (j^aXj^iXj ^ ~ b^ XjXk^J 

-2£ a J dd R fi (daXjdiXj + ^ijkdaXjXk) , (72) 


where we have used Greens theorem, and r is the unit vector in the radial 
direction. 

At this point we may substitute the asymptotic form of Xi- Ignoring the 
2 dependence 

Xi -> Xo^y and djXi (^~ ~ (73) 

Inserting this into the equation above we find that the linear term identically 
vanishes as required. 

We now proceed to calculate the second order contribution to the energy 
density 

1 d 2 V 

6 2 E = djdxidjdxi + r]£ijkSxidjSxk + 2 dy dy SxiS X j 

( 1 g’ 2 y 

= dj (SxidjSxi) + Sxi ( -d 2 Sxi + r]£ijkdj5xk + 2 d\ i)\ * Xj 

(74) 



Inserting 6\i = -£ad a Xi 


5 2 E 


fia£,bdi (daXjdidbXj) 


4“ Ca £ b da X i &b 


(- Q2 Xi + V£ijkd jX k + 


1 dV\ 

2 dxi) 


£a£,bdi (daXjdidbXj) j 


(75) 


where we have used the equation of motion in the last line. 

We may now integrate over the perpendicular coordinates and use Green’s 
theorem: 


5 2 E £ a £, b / d 2 x ± di ( daXjdidbXj) = / dBR r\d a Xjd l d b x 3 (76) 


Finally, we insert the asymptotic form of x.i hi (73) and find 


^Xo, 


5'E = 


(77) 
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This particular value of m 2 was determined from our particular choice of 
boundary contour. Depending on the boundary conditions, this term may 
be modified and may even be tensorial depending on the symmetries of the 
contour choice. However, by dimensional analysis we are guaranteed to have 

( 78 ) 
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